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Abstract 



A body dissipates energy when it freely rotates about any axis different from prin- 
cipal. This entails relaxation, i.e., decrease of the rotational energy, with the angular 
momentum preserved. The spin about the major-inertia axis corresponds to the minimal 
kinetic energy, for a fixed angular momentum. Thence one may expect comets and aster- 
oids (as well as spacecraft or cosmic-dust granules) stay in this, so-called principal, state 
of rotation, unless they are forced out of this state by a collision, or a tidal interaction, or 
cometary jetting, or by whatever other reason. As is well known, comet P/Halley, asteroid 
4179 Toutatis, and some other small bodies exhibit very complex rotational motions at- 
tributed to these objects being in non-principal states of spin. Most probably, the asteroid 
and cometary wobble is quite a generic phenomenon. The theory of wobble with internal 
dissipation has not been fully developed as yet. In this article we demonstrate that in 
some spin states the effectiveness of the inelastic-dissipation process is several orders of 
magnitude higher than believed previously, and can be measured, by the presently avail- 
able observational instruments, within approximately a year span. We also show that in 
some other spin states both the precession and precession-relaxation processes slow down 
considerably. (We call it near-separatrix lingering effect.) Such spin states may evolve so 
slowly that they can mimic the principal-rotation state. 
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I. INTRODUCTION 



It was a surprise for mission experts when, in 1958, the Explorer satellite changed its rotation axis. 
The satellite, a very prolate body with four deformable antennas on it, was planned to spin about its 
least-inertia axis, but for some reason refused to do so. Later the reason was understood: on general 
grounds, the body should end up in the spin state that minimises the kinetic rotational energy, for a fixed 
angular momentum. The rotation state about the maximal-inertia axis is the one minimising the energy, 
whereas spin about the least-inertia axis corresponds to the maximal energy. As a result, the body must 
get rid of the excessive energy and change the spin axis. This explains the vicissitudes of the Explorer 
mission. 

Similarly to spacecraft, a comet or an asteroid in a non-principal rotation mode will dissipate energy 
and will, accordingly, return to the stable spin (Black et al. 1999, Efroimsky & Lazarian 1999). Never- 
theless, several objects were found in excited states of rotation. These are, for example, comet P/Halley 
(Sagdeev et al. 1989; Peale & Lissauer 1989; Peale 1991), comet 46P/Wirtanen (Samarasinha, Mueller 
& Belton 1996; Rickman & Jorda 1998), comet 29P/Schwachmann-Wachmann 1 (Meech et al 1993), 
asteroid 1620 Geographos (Prokof'eva et al. 1997; Prokof'eva et al. 1996), and asteroid 4179 Toutatis 
(Ostro et al. 1993, Harris 1994, Ostro et al. 1995, Hudson and Ostro 1995, Scheeres et al. 1998, Ostro 
et al. 1999). 

The dynamics of a freely rotating body is determined, on the one hand, by the initial conditions of 
the object's formation and by the external factors forcing the body out of its principal spin state. On the 
other hand, it is influenced by the internal dissipation of the excessive kinetic energy associated with 
wobble. Two mechanisms of internal dissipation are known. The so-called Barnett dissipation, caused 
by the periodic remagnetisation, is relevant only in the case of cosmic-dust-granule alignment (Lazarian 
& Draine 1997). The other mechanism, called inelastic relaxation, is also relevant for mesoscopic 
grains, and plays a primary role in the case of macroscopic bodies 

Inelastic relaxation results from the alternating stresses that are generated inside a wobbling body 
by the transversal and centripetal acceleration of its parts. The stresses deform the body, and inelastic 
effects cause energy dissipation. 

The external factors capable of driving a rotator into an excited state are impacts and tidal interac- 
tions, the latter being of a special relevance for planet-crossers. In the case of comets, wobble is largely 
caused by jetting. Even gradual outgassing may contribute to the effect because a spinning body will 
start tumbling if it changes its principal axes through a partial loss of its mass or through some redis- 
tribution thereof. Sometimes the entire asteroid or comet may be a wobbling fragment of a progenitor 
disrupted by a collision (Asphaug & Scheeres 1999, Giblin & Farinella 1997, Giblin et al. 1998) or by 
tidal forces. All these factors that excite rotators compete with the inelastic dissipation that always tends 
to return the rotator back to the fold. 

Study of rotation of small bodies may provide much information about their recent history and 
internal structure. However, theoretical interpretation of the observational data will become possible 
only after we understand quantitatively how inelastic dissipation affects rotation. 

Evidently, the kinetic energy of rotation will decrease at a rate equal to that of energy losses in the 
material. Thus one should first calculate the elastic energy stored in a tumbling body, and then calculate 
the energy-dissipation rate, using the material quality factor Q. This empirical factor is introduced for 
a phenomenological description of the overall effect of the various attenuation mechanisms (Nowick & 
Berry 1972; Burns 1986, 1977; Knopoff 1963; Goldreich & Soter 1965). A comprehensive discussion 
of the Q-factor of asteroids and of its frequency- and temperature-dependence is presented in Efroimsky 
& Lazarian (2000). 
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A pioneer attempt to study inelastic relaxation in a wobbling asteroid was made by PrendergastQ 
back in 1958. Prendergast (1958) pointed out that the precession-caused acceleration must create fields 
of stress and strain over the body volume. He did not notice the generation of the higher-than-the-second 
harmonics, but he did point out the presence of the second harmonic: he noticed that precession with 
rate uj produces stresses of frequency 2 uj along with those of uj . Since Prendergast's calculation 
was wrong in several respects (for example, he missed the term fl x r ), it gave him no chance to 
correctly evaluate the role of nonlinearity, i.e., to estimate the relative contribution of the harmonics to 
the entire effect, a contribution that is sometimes of the leading order. Nonetheless, Prendergast should 
be credited for being the first to notice the essentially nonlinear character of the inelastic relaxation. 
Another publication where the emergence of the second harmonic was pointed out was by Peale (1973) 
who addressed inelastic relaxation in the case of nearly spherical bodies. This key observation made 
by Prendergast and Peale went unnoticed by colleagues and was forgotten, even though their papers 
were once in a while mentioned in the references. Later studies undertaken by Burns & Safronov 
(1973), for asteroids, and by Purcell (1979), for cosmic-dust granules, treated the issue from scratch 
and fully ignored not only the higher harmonics but even the second mode 2 uj . This led them to 
a several-order-of-magnitude underestimation of the effectiveness of the process, because the leading 
effect comes often from the second mode and sometimes from higher modes. One more subtlety, missed 
by everyone who ever approached this problem, was that, amazingly, the harmonics uj n = nuj\ are not 
necessarily multiples of the precession rate uj . We shall demonstrate that in fact a body precessing at 
rate uj experiences a superposition of stresses alternating at frequencies uj n = nu%. Here the "base 
frequency" uj\ can either coincide with the precession rate uj or be lower than it, dependent upon 
the symmetry of the top and upon its rotation state. (For example, uj 1 coincides with uj in the case of 
symmetrical oblate rotator. In this special case, uj\ itself and uj 2 = 2uj x = 2uj are the only emerging 
harmonics. But in the general case of a triaxial top all the other harmonics will show themselves.) 

Another oversight present in all the afore-quoted studies was their mishandling of the boundary 
conditions. In Purcell's article, where the body was modelled by an oblate rectangular prism, the normal 
stresses had their maximal values on the free surfaces and vanished in the centre of the body (instead of 
being maximal in the centre and vanishing on the surfaces). In Burns & Safronov (1973) the boundary 
conditions, at first glance, were not touched upon at all. In fact, they were addressed tacitly when the 
authors tried to decompose the pattern of deformations into bending and bulge flexing. An assumption 
adopted in Burns & Safronov (1973), that "the centrifugal bulge and its associated strains wobble back 
and forth relative to the body as the rotation axis uj moves through the body during a wobble period," 
lead the authors to a misconclusion that the "motion of the bulge through the (nutation) angle a 
produces strain energy" and to a calculation based thereon. In reality, however, the bulge appearance is 
but an iceberg tip, in that an overwhelming part of the inelastic dissipation process is taking place not 
near the surface but in the depth of the body, deep beneath the bulge. This follows from the fact that the 
stress and strain are small in the shallow regions and increase in the depth, if the boundary conditions 
are improved. (Remarkably, in the Peale, Cassen and Reynolds (1979) discussion of tidal dissipation in 
Io, maximum dissipation was also found to occur in the centre of the initially solid Io.) 



We wish to thank Vladislav Sidorenko for drawing our attention to Prendergast's article, and for reminding us 
that much that appears new is well-forgotten past. 

2 Term suggested by William Newman. 
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II. NOTATIONS AND ASSUMPTIONS 



We shall use two coordinate systems. The body frame will naturally be represented by the three 
principal axes of inertia: 1, 2, and 3, with coordinates x, y, z, and unit vectors d, e 2 , e 3 . The second 
(inertial) frame (X, Y, Z), with basis vectors ex, ey, ez, may be chosen with its Z axis aimed along 
the body's angular-momentum vector J and with its origin coinciding with that of the body frame (i.e., 
with the centre of mass). The inertial-frame coordinates will be denoted by the same capital letters: X, 
Y, and Z as the axes. 

The angular momentum J of a freely -preces sing body is conserved in the inertial frame. The body- 
frame-related components of the inertial angular velocity fl will be called ^1,2,3 , while letter lu will 
be reserved for the precession rate. 

A body-frame-based observer will view both the inertial angular velocity fi and the angular mo- 
mentum J nutating around the principal axis 3 at rate lu . An inertial observer, though, will argue that 
it is rather axis 3 and angular velocity O that are wobbling about J. As is well known, the precession 
rate of O about J in the inertial frame is different from the precession rate lu of Q about axis 3 
in the body frame. (See Section IV.) We would emphasize that the precession rate of our interest is the 
one in the body frame. 

Free rotation of a body is described by Euler's equations 

j t P,) = (Ij - I k ) Q,Q k , (2.1) 
(ijk) being a cyclic transposition of (123), and the principal moments of inertia ranging as 

h < h < h ■ (2.2) 

Since no uniform concensus on notations exists in the literature, the following table may simplify read- 
ing: 
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TABLES 





Principal moments 
of inertia 


Components of the 
angular 

velocity in the body 
frame 


Rate of angular- 
velocity wobble, 
in the body frame 


(Purcell 1979), 
(Lazarian & 
Efroimsky 1999), 
(Efroimsky 2000), 
(Efroimsky & 
Lazarian 2000), 
present article 


h > h > h 






(Synge & Griffiths 
1959) 


A>B>C 


^1 , ^2 , ^3 


P 


(Black et al. 1999) 


C>B>A 




V 
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In the body frame, the period of angular- velocity precession about the principal axis 3 is: r = 
2 7r juj . Evidently, 



aM « r- 1 , h/h « r^e , (2.3) 

e being a typical value of the relative strain that is several orders less than unity. These estimates lead 
to the inequality Ii Vti <C Ii&i , thereby justifying the commonly used approximation to Euler's 
equations: 

h tk = (Ij - h) V k . (2.4) 

Thus it turns out that in our treatment the same phenomenon is neglected in one context and accounted 
for in another: on the one hand, the very process of the inelastic dissipation stems from the precession- 
inflicted small deformations; on the other hand, we neglect these deformations in order to write down 
fl2.4p. This approximation (also discussed in Lambeck 1988) may be called adiabatic, and it remains 
acceptable insofar as the relaxation is slow against rotation and precession. To cast the adiabatic ap- 
proximation into its exact form, one should first come up with a measure of the relaxation rate. Clearly, 
this should be the time derivative of the angle 9 made by the major-inertia axis 3 and the angular 
momentum J . The axis aligns towards J , so 9 must eventually decrease. Be mindful, though, 
that even in the absence of dissipation, 9 does evolve in time, as can be shown from the equations 
of motion. Fortunately, this evolution is periodic, so one may deal with a time derivative of the angle 
averaged over the precession period. In practice, it turns out to be more convenient to deal with the 
squared sine of 9 (Efroimsky 2000) and to write the adiabaticity assertion as: 



d (sin 2 

uj being the precession rate and < ... > being the average over the precession period. The case of an 
oblate symmetrical topQ is exceptional, in that 9 remains, when dissipation is neglected, constant over 
a precession cycle. No averaging is needed, and the adiabaticity condition simplifies to: 




(2.6) 

(oblate) 

We would emphasise once again that the distinction between the oblate and triaxial cases, distinction 
resulting in the different forms of the adiabaticity condition, stems from the difference in the evolution of 
9 in the weak-dissipation limit. The equations of motion of an oblate rotator show that, in the said limit, 
9 stays virtually unchanged through a precession cycle (see section IV below). So the slow decrease 
of 9 , accumulated over many periods, becomes an adequate measure for the relaxation rate. The rate 
remains slow, compared to the rotation and precession, insofar as ( |2.6D holds. In the general case of a 
triaxial top the equations of motion show that, even in the absence of dissipation, angle 9 periodically 
evolves, though its average over a cycle stays unchanged (virtually unchanged, when dissipation is 
present but weak)^. In this case we should measure the relaxation rate by the accumulated, over many 



3 Hereafter oblate symmetry will imply not a geometrical symmetry but only the so-called dynamical symmetry: 

h = h. 

4 See formulae (Al) - (A4) in the Appendix to Efroimsky 2000. 
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cycles, change in the average of 9 (or of sin 2 9 ). Then our assumption about the relaxation being slow 
yields (ET 



The above conditions (2.5) - (|2.6[) foreshadow the applicability domain of our further analysis. For 



example, of the two quantities, 

J 2 fi 2 + J 2 fi 2 + Jffi 2 = J 2 , (2.7) 

I\ Vt\ + J 2 ^2 + -^3^3 = 2 Tfcj n , (2.8) 

only the former will conserve exactly, while the latter will remain virtually unchanged through one cycle 
and will be gradually changing through many cycles (just like (sm 2 6^ ). 

IE. THE STRATEGY 

As mentioned above, in the case of an oblate body, when the moments of inertia relate as J 3 > 
I 2 = 1 1 , the angle 9 between axis 3 and J remains adiabatically unchanged over the precession 
cycle. Hence in this case we shall be interested in 9, the rate of the maximum-inertia axis' approach to 
the direction of J. In the general case of a triaxial rotator, angle 9 evolves through the cycle, but its 
evolution is almost periodic and, thus, its average over the cycle remains virtually constant. Practically, 
it will turn out to be more convenient to use the average of its squared sine. In this case, the alignment 
rate will be characterised by the time derivative of < sin 2 9 > . Evidently, 

d < sin 2 9 > d < sin 2 9 > dT kin 
dt dT kin dt 

while for an oblate case, when 9 remains virtually unchanged over a cycle, one would simply write: 

fd9\ = ( de \ dT km 

The derivative d < sin 2 9 > /dT kin appearing in (|3TT|), as well as (d 6 / dT kin ) r Mate \ appearing in 



(|3.2[), can be calculated from Q2.7[), Q2.8D and the equations of motion. These derivatives indicate how 
the losses of the rotational energy affect the value of < sin 2 9 > (or simply of 9, in the oblate case). 
The kinetic-energy decrease is caused by the inelastic dissipation, 

f km =<W> , (3.3) 

W being the elastic energy of the alternating stresses, and < W > being its average over a precession 
cycle. This averaging is justified within our adiabatic approach. So we shall eventually deal with the 
following formulae for the alignment rate: 



d < sin 9 > d < sin 9 > d <W > 



dt dT kin dt 



(3.4) 



in the general case, and 



'd9\ ( d9 \ d <W > 



for an oblate rotator. 
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Now we are prepared to set out the strategy of our further work. While calculation of 
d( sin 2 9 ) /dTkin and (d6 / 'dTkin) ob i ate is an easy exerciseQ our main goal will be to find the dis- 
sipation rate d < W > jdt . This quantity will consist of inputs from the dissipation rates at all 
the frequencies involved in the process, i.e., from the harmonics at which stresses oscillate in a body 
precessing at a given rate to . The stress is a tensorial extension of the notion of a pressure or force. 
Stresses naturally emerge in a spinning body due to the centripetal and transversal accelerations of its 
parts. Due to the precession, these stresses contain time-dependent components. If we find a solution to 
the boundary-value problem for alternating stresses, it will enable us to write down explicitly the time- 
dependent part of the elastic energy stored in the wobbling body, and to separate contributions from 
different harmonics: 

< W > = ]T < W(u n ) > . (3.6) 

n 

W(u n ) being the elastic energy of stresses alternating at frequency to n . One should know each con- 
tribution W(to n ), for these will determine the dissipation rate at the appropriate frequency, through the 
frequency-dependent empirical quality factors. The knowledge of these factors, along with the averages 
< W(to n ) > , will enable us to find the dissipation rates at each harmonic. Sum of those will give the 
entire dissipation rate due to the alternating stresses emerging in a precessing body. 



IV. INELASTIC DISSIPATION 

Equation (|3.6[ ) implements the most important observation upon which all our study rests: gener- 
ation of harmonics in the stresses inside a precessing rigid body. The harmonics emerge because the 
acceleration of a point inside a precessing body contains centrifugal terms that are quadratic in the angu- 
lar velocity f2. In the simpliest case of a symmetrical oblate body, for example, the body-frame-related 
components of the angular velocity are given in terms of sin tot and cos tot (see formulae (p.2[) from 



Section V). Evidently, squaring of CI will yield terms both with sin cut or cos cut and with sin2u;t or 
cos 2ut . The stresses produced by this acceleration will, too, contain terms with frequency tot as well 
as those with the harmonic 2uot. In the further sections we shall explain that a triaxial body precessing at 
rate to is subject, in distinction from a symmetrical oblate body, to a superposition of stresses oscillating 
at frequencies to n = nui, the "base frequency" to x being lower than the precession rate to. The basic 
idea is that in the general, non-oblate case, the time dependence of the acceleration and stresses will 
be expressed not by trigonometric but by elliptic functions whose expansions over the trigonometric 
functions will generate an infinite number of harmonics. 

The total dissipation rate will be a sum of the particular rates (Stacey 1992) to be calculated em- 
pirically. The empirical description of attenuation is based on the quality factor Q(to) and on the 
assumption of attenuation rates at different harmonics being independent from one another: 

w = E *K) = - £ ^ = - 2 E — > < 4 -d 

Q(to) being the quality factor of the material, and Wo(u n ) and < W(to n ) > being the maximal and 
the average values of the appropriate-to-c<j n fraction of elastic energy stored in the body. This expression 
will become more general if we put the quality factor under the integral, implying its possible coordinate 
dependence]]: 



5< 



See formula ( |5.18| ) below and also formulae (A12 - A 13) in Efroimsky 2000. 



6 In strongly inhomogeneous precessing bodies the attenuation may depend on location. 
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The above assumption of attenuation rates at different harmonics being mutually independent is justified 
by the extreme smallness of strains (typically, much less than 1CT 6 ) and by the frequencies being 
extremely low (10 -5 — 1CT 3 Hz). One, thus, may say that the problem is highly nonlinear, in that we 
shall take into account the higher harmonics in the expression for stresses. At the same time, the problem 
remains linear in the sense that we shall neglect any nonlinearity stemming from the material properties 
(in other words, we shall assume that the strains are linear function of stresses). We would emphasize, 
though, that the nonlinearity is most essential, i.e., that the harmonics to n come to life unavoidably: no 
matter what the properties of the material are, the harmonics do emerge in the expressions for stresses. 
Moreover, as we shall see, the harmonics interfere with one another due to W being quadratic in 
stresses. Generally, all the infinite amount of multiples of Ui will emerge. The oblate case, where only 
lui and 2ui show themselves, is an exception. Another exception is the narrow-cone precession of a 
triaxial rotator studied in Efroimsky (2000): in the narrow-cone case, only the first and second modes 
are relevant (and uj\ ps uj). 

Often the overall dissipation rate, and therefore the relaxation rate is determined mostly by har- 
monics rather than by the principal frequency. This fact was discovered only recently (Efroimsky & 
Lazarian 2000, Efroimsky 2000, Lazarian & Efroimsky 1999), and it led to a considerable re-evaluation 
of the effectiveness of the inelastic-dissipation mechanism. In some of the preceding publications, its 
effectiveness had been underestimated by several orders of magnitude, and the main reason for this 
underestimation was neglection of the second and higher harmonics. As for the choice of values of 
the quality factor Q , Prendergast (1958) and Burns & Safronov (1973) borrowed the terrestial seis- 
mological data for Q. In Efroimsky & Lazarian (2000), we argue that these data are inapplicable to 
asteroids. 

To calculate the afore mentioned average energies < W(u n ) > , we use such entities as stress and 
strain. As already mentioned above, the stress is a tensorial generalisation of the notion of pressure. The 
strain tensor is analogous to the stretching of a spring (rendered in dimensionless fashion by relating the 
displacement to the base length). Each tensor component of the stress consists of two inputs, elastic and 
plastic. The former is related to the strain through the elasticity constants of the material; the latter is 
related to the time-derivative of the strain, through the viscosity coefficients. As our analysis is aimed at 
extremely small deformations of cold bodies, the viscosity may well be neglected, and the stress tensor 
will be approximated, to a high accuracy, by its elastic part. Thence, according to Landau & Lifshitz 
(1976), the components of the elastic stress tensor are interconnected with those of the strain tensor 
€y like: 



e *i = S H 1^17 + a v ~ ~ 6 H Tra \— (4.3) 



Tra ( 1 \ 1 

9K \ v 3 13 J 2p 

p and K being the adiabatic shear and bulk moduli, and Tr standing for the trace of a tensor. 

To simplify the derivation of the stress tensor, the body will be modelled by a rectangular prism of 
dimensions 2ax26x2c where a > b > c. The tensor is symmetrical and is defined by 

diOij = paj , (4.4) 

dj being the time-dependent parts of the acceleration components, and paj being the time-dependent 
parts of the components of the force acting on a unit volume[|. Besides, the tensor must obey the 



Needless to say, these acceleration components aj are not to be mixed with a which is the longest dimension 
of the prism. 
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boundary conditions: its product by normal unit vector, a^rij , must vanish on the boundaries of the 
body (this condition was not fulfilled in Purcell (1979)). 

Solution to the boundary-value problem provides such a distribution of the stresses and strains over 
the body volume that an overwhelming share of dissipation is taking place not near the surface but in the 
depth of the body. For this reason, the prism model gives a good approximation to realistic bodies. Still, 
in further studies it will be good to generalise our solution to ellipsoidal shapes. Such a generalisation 
seems to be quite achievable, judging by the recent progress on the appropriate boundary- value problem 
(Denisov & Novikov 1987). 

Equation ( ft~4] ) has a simple scalar analogue^ Consider a non-rotating homogeneous liquid planet of 
radius R and density p. Let g(r) and P(r) be the free-fall acceleration and the self-gravitational 
pressure at the distance r < R from the centre. (Evidently, g(r) = (4/3) n G pr .) Then the analogue 



to (4.4) will read: 



P9{r) = -?^ , (4.5) 

the expression pg(r) standing for the gravity force acting upon a unit volume, and the boundary 
condition being P(R) = 0. Solving equation ( |Q1 ) reveals that the pressure has a maximum at the 
centre of the planet, although the force is greatest at the surface. Evidently, the maximal deformations 
(strains) also will be experienced by the material near the centre of the planet. 

In our case, the acceleration a of a point inside the precessing body will be given not by the free-fall 
acceleration g(r) but will be a sum of the centripetal and transversal accelerations: Q x (H x r) + SI x 
r , the Coriolis term being negligibly small. Thereby, the absolute value of a will be proportional to 
that of r , much like in the above example. In distinction from the example, though, the acceleration of a 
point inside a wobbling top will have both a constant and a periodic component, the latter emerging due 
to the precession. For example, in the case of a symmetrical oblate rotator, the precessing components of 
the angular velocity ft will be proportional to sin tut and cos cut , whence the transversal acceleration 
will contain frequency u while the centripetal one will contain 2cu. The stresses obtained through (|4.4|) 



will oscillate at the same frequencies, and so will the strains. As we already mentioned, in the case of 
a non-symmetrical top an infinite amount of harmonics will emerge, though these will be obertones not 
of the precession rate cu but of some different "base frequency" Ui that is less than u. 

Here follows the expression for the (averaged over a precession period) elastic energy stored in a 
unit volume of the body: 



d <W > 1 1 f / 2 fi 1 



- < €ij > 



dV 2 J J 4/i IA9 K 

1 



< (TV a) 2 > + <aijaij>\ 



4/i I 1 + v 



T < (TraY > + < a 2 x > + < a 2 y > + < a 2 zz > +2 < a 2 xy + a 2 z + a 2 zx > (4.6) 



where 2p,/(9K) — 1/3 = — vj{\ + v) ~ —1/5, v being Poisson's ratio (for most solids v ps 1/4). 
Naturally^, the total averaged elastic energy is given by the integral over the body's volume: 

<W>= i J dVa i:j eij , (4.7) 



8 This very simple and, nonetheless, so illustrative example was kindly offered to me by William Newman. 

9 Very naturally indeed, because, for example, a xx e xx dV = (a xx dy dz)(e xx dx) is a product of the x-directed 
pressure upon the x-directed elongation of the elementary volume dV. 
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and it must be expanded into the sum ( |3T6| ) of inputs from oscillations of stresses at different frequencies. 
Each term (W(u n )) emerging in that sum will then be plugged into the expression (|4.ip, together with 
the value of Q appropriate to the overtone cu n . 



V. A SPECIAL CASE: PRECESSION OF AN OBLATE BODY. 



An oblate body has moments of inertia that relate as: 



h > h 



h = I 



(5.1) 



We shall be interested in 9, the rate of the maximum-inertia axis' approach to the direction of angular 
momentum J. To achieve this goal, we shall have to know the rate of energy losses caused by the peri- 
odic deformation. To calculate this deformation, it will be necessary to find the acceleration experienced 
by a particle located inside the body at a point (x, y, z). Note that we address the inertial acceleration, 
i.e., the one with respect to the inertial frame (X, Y, Z), but we express it in terms of coordinates x, y 
and z of the body frame (1, 2, 3) because eventually we shall have to compute the elastic energy stored 
in the entire body (through integration of the elastic-energy density over the body volume). 

The fast motions (revolution and precession) obey, in the adiabatical approximation, the simplified 
Euler equations (|2.4[). Their solution, with neglect of the slow relaxation, looks (Fowles and Cassiday 



1986, Landau and Lifshitz 1976), in the oblate case (5.1) 



Six 



Q i cos tut , 



Q i sin uot , 



const 



(5.2) 



where 



sm a 



fi 3 = ft 



cos a 



(5.3) 



a being the angle made by the major-inertia axis 3 with f2 . Expressions ( p72\ ) show that in the body 
frame the angular velocity f2 describes a circular cone about the principal axis 3 at a constant rate 



u = (h - l)fi 3 , h = hj I 



(5.4) 



So angle a remains virtually unchanged through a cycle (though in the presence of dissipation it still 
may change gradually over many cycles). The precession rate uj is of the same order as |f2|, except in 
the case of h — > 1 or in a very special case of Ct and J being orthogonal or almost orthogonal to the 
maximal-inertia axis 3. Hence one may call not only the rotation but also the precession "fast motions" 
(implying that the relaxation process is a slow one). Now, let 9 be the angle between the principal axis 
3 and the angular-momentum J , so that J 3 = J cos 9 and 



wherefrom 



^3 = — - = — cos 9 



cu = (h — 1) — cos6> 



(5.5) 



(5.6) 



Since, for an oblate object, 



ii fii ex + 7 2 Q 2 e 2 + h ^3 e 3 = / (^i ex + Q 2 e 2 ) + h ^3 e 3 



(5.7) 



the quantity Vt^ = JVt\ + fi 2 ^ s connected with the absolute value of J like: 
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Q i 



J 



sin 6 1 



J 



h sin 6 1 , h = 



(5.8) 



It ensues from ( |5.3D that fi_|_/n 3 
ft, tan 9 . Hence, 



tan a . On the other hand, ( |5.5D and ( |5.8| ) entail: 



tana = h tan# 



(5.9) 



We see that angle 9 is almost constant too (though it gradually changes through many cycles). We also 
see from ( |5T7| ) that in the body frame the angular-momentum vector J describes a circular cone about 
axis 3 with the same rate w as O . An inertial observer, though, will insist that it is rather axis 3, as 



well as the angular velocity Q , that is describing circular cones around J. It follows trivially from ( |5.4| ) 
and fl577p that 



I n 



I uj e 3 



(5.10) 



whence it is obvious that, in the inertial frame, both ft and axis 3 are precessing about J at rate J/ 1. 
(The angular velocity of this precession is e 3 = x e 3 = (J// — uj e 3 ) x e 3 = (3/1) x e 3 .) 
Interestingly, the rate uj = (h — 1) Q 3 , at which fl and J are precessing about axis 3 in the body 
frame, differs considerably from the rate J/ 1 at which ft and axis 3 are precessing around J in the 
inertial frame. (In the case of the Earth, J/ 1 « 400 uj because h is close to unity.) Remarkably, 
the inertial-frame-related precession rate is energy-independent and, thus, stays unchanged through the 
relaxation process. This is not the case for the body-frame-related rate uj which, according to (|5.6|), 
gradually changes because so does 9. 

As is explained above, we shall be interested in the body-frame-related components ^1,2,3 precessing 
at rate uj about the principal axis 3. Acceleration of an arbitrary point of the body can be expressed in 
terms of these components through formula 



fi x 



2 ft x v' + fi x (fi x r') 



(5.11) 



where r, v, a are the position, velocity and acceleration in the inertial frame, and r', v' and a' 
are those in the body frame. Here r = r' . Mind though that v' and a' do not vanish in the 
body frame. They may be neglected on the same grounds as term IiVti in (|2TT|): precession of a body 
of dimensions ~ I , with period r , leads to deformation-inflicted velocities v' ~ el/r and 



accelerations 
very small e 



el/r' 



being the typical order of strains arising in the material. Clearly, for 



quantities v' and a' are much less than the velocities and accelerations of the body as 



a whole (that are about l/r and l/r 2 , correspondingly). Neglecting these, we get, from ( |5.1 1[ ) and 



<K3b, for the acceleration at point (x, y, z) 



+ e 2 I - Q 2 ± x sin 2ujt — — Q 2 ± y cos 2ujt + z Q±_ fi 3 h sin ut X + 



+ e 3 ^3 (2 — h) (x cos ujt + y sin ujt ) } 



(5.12) 



Plugging this into Q4.4D, with the proper boundary conditions imposed, yields, for an oblate prism of 
dimensions 2a x 2a x 2c , a > c : 



[X 



cos 2ut 



a 



yy 



{y 1 



cos 2ujt , a z 







(5.13) 
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a 



■HI 



n 2 ± 



' 2 i 2 

x + y 



2a ) sin 2ut 



(5.14) 



C7, 



/i (2 2 - c 2 ) + (2 - h) {x 2 



cos cut 



(5.15) 



Ol tt 3 h {z 2 - c 2 ) + (2 - h) (y 2 



sin out 



(5.16) 



In ( |5.12D - ( |5.16| ) we kept only time-dependent parts, because time-independent parts of the acceleration, 
stresses and strains are irrelevant in the context of dissipation. A detailed derivation of Q5.12D - ( p~.16[ ) is 
presented in (Lazarian & Efroimsky 1999). 

Formulae ( |5.13[ ) - ( |5.16| ) implement the polynomial approximation to the stress tensor. This approx- 
imation keeps the symmetry and obeys Q4.4| ) with ( |5.12| ) plugged into it. The boundary condition are 
satisfied exactly for the diagonal components and only approximately for the off-diagonal components. 
The approximation considerably simplifies calculations and entails only minor errors in the numerical 
factors in ( frZip . 

The second overtone emerges, along with the principal frequency to , in the expressions for stresses 
since the centripetal part of the acceleration is quadratic in Q . The kinetic energy of an oblate spinning 
body reads, according to (|2.8[), Q5.3|), and (|5.9[): 



T, 



kin 



in 2 



wherefrom 



dTi 



kin 



J 2 

f (h 



1 
2 



- sin 2 6 + — cos 2 6 



1) sin 9 cos 9 



uj J sin 9 



The latter expression, together with ( [3.5D and (rkq), leads to 



d9 
It 



d9 



kin 



(It 



(oj J sin 



W , 



J 



(5.17) 



(5.18) 



(5.19) 



where 

W = W {lv) + W {2uj) = uj 



u) 



+ 2uj 



W, 



(2w) 



2oj 



Q(2u J ) 



— {< W {u) > + 2 < W i2u}) >} , (5.20) 



the quality factor assumed to depend upon the frequency very weakly^. In the above formula, Wq and 
are amplitudes of elastic energies corresponding to the principal mode and the second harmonic. 
Quantities < W w >= Wq/2 and < W 2ul >= W^/2 are the appropriate averages. Substitution of 



( J5 . 1 3[ ) - ( |5.16| ) into (|4.6|), with further integration over the volume and plugging the result into (|3.5D, will 
give us the final expression for the alignment rate: 



d9/dt 



sin 3 ^ 

2 4 



63 (c/a) 4 cot 2 9 + 20 a 2 n 3 p 



[1 + (c/ 



a 



(214 



H Q 



(5.21) 



10 The uj -dependence of Q should be taken into account within frequency spans of several orders, but is 
irrelevant for frequencies differing by a factor of two. 
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where 



= t (5.22) 

^3 

is a typical angular velocity. Deriving ( |5.21| ), we took into account that, for an oblate 2a x 2a x 2c 
prism (where a > c ), the moment of inertia J 3 and the parameter /i read: 

^s = ^-pac , h = — = . (5.23) 

3 7 1+ {c/a) 2 



Details of derivation of ( |5.21[ ) are presented in (Lazarian & Efroimsky 1999)0- 



Formula Q5.21D shows that the major-inertia axis slows down its alignment at small residual angles. 
For 9 — > 0, the derivative 9 becomes proportional to 9, and thus, 9 decreases exponentially slowly: 
9 = A exp(— £t), where A and ( are some positive numbers^. This feature, "exponentially slow finish", 
(which was also mentioned, with regard to the Chandler wobble, in Peale (1973), formula (55)) is natural 
for a relaxation process, and does not lead to an infinite relaxation time if one takes into account the 
finite resolution of the equipment. Below we shall discuss this topic at length. 

Another feature one might expect of ( p.21[ ) would be a "slow start": it would be good if d9/dt 
could vanish for 9 — ► ir/2 . If this were so, it would mean that at 9 = n/2 (i.e., when the major- 
inertia axis is exactly perpendicular to the angular-momentum vector) the body "hesitates" whether to 
start aligning its maximal-inertia axis along or opposite to the angular momentum, and the preferred 
direction is eventually determined by some stochastic influence from the outside, like (say) a collision 
with a small meteorite. This behaviour is the simplest example of the famous spontaneous symmetry 
breaking, and in this setting it is desirable simply for symmetry reasons: 9 = n/2 must be a position 



of an unstable equilibrium^. Contrary to these expectations, though, Q5.21| ) leaves d9/dt nonvanishing 



for 9 — ► 7r/ 2 , bringing an illusion that the major axis leaves the position 9 = n/2 at a finite rate. This 
failure of our formula ( |5.21| ) comes from the inapplicability of our analysis in the closemost vicinity of 
9 = 7r/2 . This vicinity simply falls out of the adiabaticity realm adumbrated by (pT6|), because u 
given by (|5.6D vanishes for 9 — > ir/2 (then one can no longer assume the relaxation to be much slower 
than the precession rate, and hence, the averaging over period becomes illegitimate). 

One more situation, that does not satisfy the adiabaticity assertion, is when uj vanishes due to (h — 
1) — > 0. This happens when cj a approaches unity. According to (|5.21[), it will appear that d9 / dt remains 



Our expression ( 5.21 ) presented here differs from the appropriate formula in Lazarian & Efroimsky (1999) by 
a factor of 2, because in Lazarian & Efroimsky (1999) we missed the coefficient 2 connecting Wq"^ with W^"'\ 

12 This resembles the behaviour of a pendulum: if the pendulum is initially given exactly the amount of kinetic 
energy sufficient for the pendulum to move up and to point upwards at the end of its motion, then formally it takes 
an infinite time for the pendulum to stand on end. 

13 Imagine a knife freely rotating about its longest dimension, and let the rotation axis be vertical. This rotation 
mode is unstable, and the knife must eventually come to rotation about its shortest dimension, the blade being in 
the horizontal plane. One cannot say, though, which of the two faces of the blade will look upward and which 
downward. This situation is also illustrated by the pendulum mentioned in the previous footnote: when put upside 
down on its end, the pendulum "hesitates" in what direction to start falling, and the choice of direction will be 
dictated by some infinitesimally weak exterior interaction (like a sound, or trembling of the pivot, or an evanescent 
flow of air). 
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nonvanishing for c/a — > 1, though on physical grounds the alignment rate must decay to zero because, 
for c = a, the body simply lacks a major- inertia axis. 

All in all, Q5.21D works when 9 is not too close to ir/2 and c/a is not too close to unity: 

J 1 - (c/a) 2 

- 9 < (h - 1) -r cos9 = - Q cos9 . (5.24) 

7 3 1 + (c/a) 2 

Knowledge of the alignment rate 9 as a function of the precession-cone half-angle 9 enables one 
not only to write down a typical relaxation time but to calculate the entire dynamics of the process. 
In particular, suppose the observer is capable of measuring the precession-cone half-angle 9 with an 



error S . This observer will then compute, by means of ( |5.21D , the time needed for the body to change 
its residual half-angle from 9 to 9 — A9 , for A 9 > 5 . This time will then be compared with 
the results of his further measurements. Below we shall show that such observations will soon become 
possible for spacecraft. 

First, let us find a typical relaxation time, i.e., a time span necessary for the major-inertia axis to 
shift considerably toward alignment with J . This time may be defined as: 

t [ 5 d ° (5 25) 

Je dO/dt ' 

#o being the initial half-angle of the precession cone (6q < tt/2), and 5 being the minimal 
experimentally-recognisable value of 9 . A finite 5 will prevent the "slow-finish" divergency. A 
particular choice of 9 and 5 will lead to an appropriate numerical factor in the final expression for 
t r . Fig.l shows that t r is not very sensitive to the choice of angle 9 , as long as this angle is not 
too small. This weak dependence upon the initial angle is natural since our approach accounts for the 
divergence at small angles ("exponentially slow finish") and ignores the "slow start". Therefore one can 
take, for a crude estimate, 

#o = vr/2 . (5.26) 

For t r it would give almost the same result as, say, 7r/3 or tt/4 . A choice of 5 
must be determined exclusively by the accuracy of the observation technique: 5 is such a min- 
imally recognizable angle that precession within a cone of half-angle 5 or less cannot be de- 
tected. Ground-based photometers measure the lightcurve-variation amplitude that is approximately 
proportional to the variation in the cross-sectional area of the wobbling body. In such sort of ex- 
periments the relative error is around 0.01. In other words, only deviations from one revolu- 
tion to the next exceeding 0.01 mag may be considered real. This corresponds to precession- 
cone half-angles 5 ~ 10° or larger (Steven Ostro, private communication). Ground-based 
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FIGURES 




FIG. 1. Precession of an oblate body: relaxation time t r as a function of 9q , where 9q is the initial value 
of the precession-cone half-angle 9 . Precession begins with the precession-cone half-angle 9 equal to 9q , and 
effectively ends when 9 reaches the minimal measurable value 5 . On this plot, it is assumed that the initial value 
9q varies from 6 = 6° through tt/2 . This choice of 5 corresponds to the current abilities of ground-based 
radars. (Spacecraft-based instruments provide S = 0.01° .) 



radars have a much sharper resolution and can grasp asteroid-shape details as fine as 10 m . This 
technique may reveal precession at half-angles of about 5 degrees. NEAR- type missions potentially may 
provide an accuracy of 0.01° (Miller et al. 1999). For a time being, we would lean to a conservative 
estimate 

5 = 6° , (5.27) 

though we hope that within the coming years this limit may be reduced by three orders due to advances 
in the spacecraft-borne instruments. 

Together, ( ]5.21[ ), ( |5.25[ ) - Q5.27D yield dependences illustrated by Fig.l. Remarkably, t r is not 



particularly sensitive to the half-sizes' ratio c/a when this ratio is between 0.5 - 0.9 (which is the case 
for realistic asteroids, comets and many spacecraft). Our formulae give: 

Hour result) « (1 - 2) -^-3 for 6 « (2 - 3) 5 = 12 - 18° ; 

P Ct i Lq 



v Q 

p a 2 VLq 



Hour result) « (3 - 4) J ~T^ f OT 9 « 7r/4 J (5.28) 



t (our result) ~ (4 - 5) for 9q ~ 7r/2 

pa i Zq 



(Mind though that, according to ( |5.24| ), 6 should not approach ir/2 too close.) To compare our results 



with a preceding study, recall that according to Burns & Safronov (1973) 
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i '- s ' ss100 ^' <5 - 29) 

The numerical factor in Burns & Safronov's formula is about 100 for objects of small oblateness, i.e., 
for comets and for many asteroids. (For objects of irregular shapes Burns and Safronov suggested a 
factor of about 20 in place of 100.) 

This numerical factor is the only difference between our formula and that of Burns & Safronov. This 
difference, however, is quite considerable: for small residual half-angles 9 , our value of the relaxation 
time is two orders smaller than that predicted by Burns & Safronov. For larger residual half-angles, 
the times differ by a factor of several dozens. We see that the effectiveness of the inelastic relaxation 
was much underestimated by our predecessors. There are three reasons for this underestimation. The 
first reason is that our calculation was based on an improved solution to the boundary-value problem 
for stresses. Expressions Q5.13D - ( |5.16D show that an overwhelming share of the deformation (and, 



therefore, of the inelastic dissipation) is taking place in the depth of the body. This is very counterin- 
tuitive, because on a heuristic level the picture of precession would look like this: a centrifugal bulge, 
with its associated strains, wobbles back and forth relative to the body as moves through the body 
during the precession period. This naive illustration would make one think that most of the dissipation 
is taking place in the shallow regions under and around the bulge. It turns out that in reality most part 
of the deformation and dissipation takes place deep beneath the bulge (much like in the simple example 
with the liquid planet, that we provided in the end of section IV). The second, most important, reason 
for our formulae giving smaller values for the relaxation time is that we have taken into account the 
second harmonic. In many rotational states this harmonic turns to be a provider of the main share of 
the entire effect. In the expression (63(c/a) 4 cot 2 9 + 20) that is a part of formula ( |5.21| ), the term 
63 (c/a) 4 cot 2 9 is due to the principal frequency, while the term 20 is due to the second harmonic^. 
For c/a belonging to the realistic interval 0.5 — 0.9 , the second harmonic contributes (after integration 
from #o through 5 ) a considerable input in the entire effect. This input will be of the leading order, 
provided the initial half-angle 9 is not too small (not smaller than about 30° ). In the case of a small 
initial half-angle, the contribution of the second mode is irrelevant. Nevertheless it is the small-angle 
case where the discrepancy between our formula and ( p.29| ) becomes maximal. The estimate ( |5.29| ) for 
the characteristic time of relaxation was obtained in Burns & Safronov (1973) simply as a reciprocal to 
their estimate for 9 ; it ignores any dependence upon the initial angle, and thus gives too long times for 
small angles. The dependence of the dissipation rate of the values of 9 is the third of the reasons for 



our results being so different from the early estimate ( |5.29D . 



Exploration of this, third, reason may give us an important handle on observation of asteroid relax- 



ation. It follows from (5.21) that a small decrease in the precession-cone half-angle, — A9 , will be 



performed during the following period of time: 

At (- A9) 2 - [1 ± (c/fl)2 J 4 1 ^ ( , , m 

{ ' 3 63 (c/a) 4 cot 2 9 + 20 sin 3 # a 2 Q 3 p ' K } 

For asteroids composed of solid silicate rock, the density may be assumed p « 2.5 x 10 3 kg/m 3 , 
while the product in the numerator should be jiQ « 1.5 x 10 13 dyne/ cm 2 = 1.5 x 10 12 Pa 
as explained in Efroimsky & Lazarian (2000). Burns & Safronov suggested a much higher value of 
3 x 10 14 dyne /cm 2 = 3 x 10 13 Pa , value acceptable within the terrestial seismology but, probably, 
inapplicable to asteroids. 



For calculational details, see Lazarian & Efroimsky (1999). 
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For asteroids composed of friable materials, Harris (1994) suggests the following values: p 
2 x 10 3 kg/m 3 and pQ m 5 x 10 12 dyne /cm 2 = 5 x 10 11 Pa . Naturally, this value is lower 
than those appropriate for solid rock (Efroimsky & Lazarian 2000), but in our opinion it is still too high 
for a friable medium. Harris borrowed the aforequoted value from preceding studies of Phobos (Yoder 
1992). Mind, though, that Phobos may consist not only of rubble: it may have a solid component in 
the centre. In this case, a purely rubble-pile asteroid may have a lower p, Q than suggested by Harris. 
Anyway, as a very conservative estimate for a rubble-pile asteroid, we shall take the value suggested by 
Harris. 

As for the geometry, let, for example, 9 = ir/3 and c/a = 0.6 . Then 

At = (-A9) 1.2 . (5.31) 

a 2 i Zg p 

If we measure time At in years, the revolution period T = 2 7r/f2 in hours, the maximal half-size 
a in kilometers, and 9 in angular degrees (A9 = A9° x 1.75 x 1CT 2 ), our formula 05.3OQ will yield: 



At {years) = (- A9°) x 1.31 x 10- 7 iiRTlp^l = 0.33 , (5.32) 

P a (km) a (km) 

where we accepted Harris' values of pQ = 5 x 10 11 Pa and p = 2 x 10 3 kg/m 3 , and the angular 
resolution of spacecraft-based devices was assumed to be as sharp as |A0| = 0.01° , according to 
Miller et al. (1999). 



VI. TRIAXIAL AND PROLATE ROTATORS 

Typically, asteroids and comets have elongated shapes, and the above formulae derived for oblate 
bodies make a very crude approximation of the wobble of a triaxial or prolate body. In the case of a 
triaxial rotator, with I 3 > I 2 > 1% , the solution to the Euler equations is expressable in terms of elliptic 
functions. According to Jacobi (1882) and Legendre (1837), it will read, for J 2 < 2 7 2 T kin , as 

Qi = 7 dn (ut, k 2 ^j , f2 2 = (3 sn (ut, k 2 ^j , Q 3 = a cn (ut, k 2 ^j , (6.1) 

while for J 2 > 2 I 2 T kin it will be: 

f2i = 7 cn (ujt, k 2 ^j , f2 2 = /3 sn (ut, k 2 ^j , = a dn (ut, k 2 ^j . (6.2) 

Here the precession rate uj and the parameters a, f3, 7 and k are some algebraic functions of 
h, 2,3, T kin and J 2 . For example, k is expressed by 



I 2 J 2 - 2hT k 



I2 — I\ 2I 3 T kin — J 2 



for (|6.1D, and by 



/ J 2 - h 2I 3 T, 



k = \It-^t ;;°^r r " , (6.4) 



I 2 J 2 - 2hT, 



kin 



for (|6.2|). In the limit of oblate symmetry (when — ► 1 ), solution ( |6.2D approaches (|5.2D, while 
the applicability region of ( |6.1| ) shrinks. Similarly, in the prolate-symmetry limit ( (I 3 — I2)/ h ) 
the applicability realm of ( |6.2D will become infinitesimally small. The easiest way of understand- 
ing this would be to consider, in the space fii , fi 2 , ^3 , the angular-momentum ellipsoid 
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J 2 = If nf + if Vt\ + If Vt\ .A trajectory described by the angular-velocity vector fl in 
the space f2i , f2 2 , fi 3 will be given by a line along which this ellipsoid intersects the kinetic-energy 
ellipsoid 2 T^ in = I\ Vt\ + J 2 fi 2 , + J 3 f2 2 , as on Fig. 2. Through the relaxation process, 
the angular-momentum ellipsoid remains unchanged, while the kinetic-energy ellipsoid evolves as the 
energy dissipates. Thus, the fast process, precession, will be illustrated by the (adiabatically) periodic 
motion of along the line of ellipsoids' intersection; the slow process, relaxation, will be illustrated 
by the gradual shift of the moving vector fl from one trajectory to another (Lamy & Burns 1972). On 
Fig. 2, we present an angular-momentum ellipsoid for an almost prolate body whose angular momenta 
relate to one another as those of asteroid 433 Eros: 1 x 3 x 3.05 (Black et al. 1999). Suppose the 
initial energy was so high that fi was moving along some trajectory close to the pole A on Fig. 2. This 
pole corresponds to rotation of the body about its minor-inertia axis. The trajectory described by f2 
about A is almost circular and remains so until fi approaches the separatrixf^]. This process will be 



described by solution (pT|). In the vicinity of separatrix, trajectories will become noticeably distorted. 




FIG. 2. The constant-angular-momentum ellipsoid, in the angular-velocity space. The lines on its surface 
are its intersections with the kinetic-energy ellipsoids corresponding to different values of the rotational energy. 
The quasi-stable pole A is the maximal-energy configuration, i.e., the state when the body spins about its mini- 
mal-inertia axis. The stable pole C symbolises the minimal-energy state, i.e., rotation about the maximal-inertia 
axis. The angular-velocity vector precesses along the constant-energy lines, and at the same time slowly shifts 
from one line to another, approaching pole C. The picture illustrates the case of an elongated body: I3 ~ I2 > I\. 
The trajectories are circular near A and (in the case of an elongated body) remain virtually circular almost up 
to the separatrix. After the separatrix is crossed (with chaotic flipovers possible), the body starts tumbling. The 
trajectories will regain a circular shape only in the closemost proximity of C. 



Crossing of the separatrix may be accompanied by stochastic flipovers[^. After the separatrix is crossed, 
librations will begin: O will be describing not an almost circular cone but an elliptic one. This process 



15 This trajectory on Fig.2 being almost circular does not necessarily mean that the precession cone of the major- 
inertia axis about J is circular or almost circular. 
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The flipovers are unavoidable if dissipation of the kinetic energy through one precession cycle is less than a 
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will be governed by solution (|6T2|). Eventually, in the closemost vicinity of pole C, the precession will 
again become almost circular. (This pole, though, will never be reached because the alignment of ft 
towards J has a vanishing rate for small residual angles.) Parameter k shows how far the tip of ft 
is from the separatrix on Fig. 2: k is zero in poles A and C, and is unity on the separatrix. It is defined 
by when ft is between pole A and the separatrix, and by ( |Q| ) when ft is between the separatrix 
and pole C. (For details see (Efroimsky 2000).) 

If in the early stage of relaxation of an almost prolate (J 3 ~ I 2 ) body the tip of vector ft is near 
pole A, then its slow departure away from A is governed by formula (9.22) in (Efroimsky 2000): 

d (sin 2 6 



dt 



4p 2 J 2 

fi Q(u) 



(h ~ h) 1 



sin 



e)) {ujS 1 



2 (sin 2 # 



111— till (| ( Mir 7, 
2 h - h h 



uj S — — ( 1 



where 



sin 2 ^)) + 2uS 2 ^4 2 -± (I - (,ur<r 



(6.5) 



'(2/ 3 T; 



kin 



J 2 ) {h ~ h) 



h h h 



JJ| (I 3 ~ h) (I 2 

h V hh 



2 ^sin 2 



(6.6) 



9 is the angle between the angular-momentum vector J and the major-inertia axis 3 ; S 0j i :2 are some 
geometrical factors (So = in the case of I 2 = h), and < ... > symbolises an average over the 
precession cycle. For < cos 2 9 > not exceeding « 1/7, this equation has an exponentially decaying 
solution. For c/a = 0.6 that solution will read: 



At 



A (9)) x 0.08 



V Q 

a 2 p 



(6.7) 



Comparing this with Q5.31D , we see that at this stage the relaxation is about 15 times faster than in the 
case of an oblate body. 

During the later stage, when ft gets close to the separatrix, all the higher harmonics will come into 
play, and our estimate will become invalid. How do the higher harmonics emerge? Plugging of ( |5.1[ ) or 
( |6.2| ) into ( |5.12| ) will give an expression for the acceleration of an arbitrary point inside the body. Due 
to ([4.4D, that expression will yield formulae for the stresses. These formulae will be similar to ( |5.13 



5.16D, but will contain elliptic functions instead of the trigonometric functions. In order to plug these 



formulae for into (|4.6|), they must first be squared and averaged over the precession cycle. For a 
rectangular prizm 2a x 2b x 2c , a direct calculation performed in (Efroimsky 2000) gives: 



a: 



.1" 



a 2 ) 2 E 1 



(6.8) 



a 



VII 



P -{S + Q) 2 p\y 2 -b 2 ) 2 Si 



(6.9) 



typical energy of an occational interaction (a tidal-force-caused perturbation, for example). 
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<£) = £ (1 - SO W - c 2 ) 2 , (6.10) 



Traf) = ?-f3 4 {(l-Q)(x 2 -a 2 ) 2 + (S + Q)(y 2 -b 2 ) + (1 - S)(z 2 - c 2 )}' E 1 , (6.11) 



<) = T {(Pl + acue^-b^ + i^-acuk 2 )^ 2 -^)}' ~ 2 , (6.12) 



4 



(<4) = J {(/*" + «7)(^ 2 - c 2 ) + (/3a; - ai )(x 2 - a 2 )} 2 ~ 3 , (6.13) 



/r f n a. ., vi 



<) = J {(«/3 + ^7)(^ 2 -c 2 ) + M-cu 7 )(y 2 -& 2 )} H 4 , (6.14) 

where Q and £ are some combinations of Ii, I 2 , h , defined by formula (2.8) in (Efroimsky 2000). 
Factors Si 2,3,4 stand for averaged powers of the elliptic functions: 



—2 



—3 



Si = (sn 2 (w, k 2 ) - < sn 2 (u, k 2 ) > 

= < sn 4 (w, k 2 ) > - < sn 2 (u, k 2 ) > 2 , (6.15) 

(^sn(u, k 2 ) cn(u, k 2 ) — sn(u, k 2 ) cn(u, 

< sn 2 (u, k 2 ) cn 2 (u, k 2 ) > - < sn(u, k 2 ) cn(u, k 2 ) > 2 , (6.16) 

( cn(u, k 2 ) dn(u, k 2 ) — ^ cn(u, k 2 ) dn(u, k 2 ) 

< cn 2 (u, k 2 ) dn 2 (u, k 2 ) > - < cn(u, k 2 ) dn(u, k 2 ) > 2 , (6.17) 



S 4 = ^ (sn(w, k 2 ) dn(u, k 2 ) — ^ sn(u, k 2 ) dn(u, k 2 ) ^ ) 
= < sn 2 (u, k 2 ) dn 2 (u, k 2 ) > — < sn(u, k 2 ) dn(u, k 2 ) > 2 , 



(6.18) 



where averaging implies: 

<...>= - ... du , (6.19) 

T JO 

t being the mutual period of sn and cn and twice the period of dn : 

r = 4K(k 2 ) = 4 (1 - k 2 sin 2 ^)" 1/2 dip . (6.20) 
Jo 
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The origin of expressions fl6.15| - 16.18| ) can be traced from formulae (8.4, 8.6 - 8.13) in (Efroimsky 



2000). For example, expression ( |5.1ip , that gives acceleration of an arbitrary point inside the body, 
contains term sn 2 (ut, k 2 ) . (Indeed, one of the components of the angular velocity is proportional to 
sn(...) , while the centripetal part of the acceleration is a quadratic form of the angular- velocity compo- 
nents.) The term sn 2 (ut, k 2 ) in the formula for acceleration yields a similar term in the expression 
for a xx . For this reason expression (8.6) in (Efroimsky 2000), that gives the time-dependent part of 
o xx , contains sn 2 {...) — < sn 2 (...)> , wherefrom ( |6.15[ ) ensues. 



Now imagine that in the formulae (6.8 - 6.14) the elliptic functions are presented by their series 



expansions over sines and cosines (Abramovitz & Stegun 1965): 

2tt 00 * q n / 2 

sn(ut, k 2 ) = —— Y 1 sm(u n t) , (6.21) 

fcA n=i L ~ 1 

2n 00 * q n l 2 

cn(ut,k 2 ) = —— Y - — - cos (u) n t) , (6.22) 
ft A n=1 i + q 

7r 2ir 00 ** q n l 2 
dn(ut, k 2 ) = — + — Y T~j — n cos (^t) , (6.23) 



where 



u n = nujn/(2K(k 2 )) , q = exp(-7rK(k' 2 )/K(k 2 )) , k' 2 = 1 - k 2 (6.24) 



and the function K{k 2 ) is the complete elliptic integral of the first kind (see Q6.20D or ( |6.29[ )). A star in 



the superscript denotes a sum over odd n's only; a double star stands for a sum over even n's. Plugging 
of ( |6.21| - |6.23| ) into ( |5T%1 - |6.14[ ) will produce, after squaring of sn, cn, dn , an infinite amount of terms 
like sin 2 (u; Tl t) and cos 2 (u> n t) , along with an infinite amount of cross terms. The latter will be removed 
after averaging over the precession period, while the former will survive for all n's and will average to 



1/2. Integration over the volume will then lead to an expression like (|3.6[), with an infinite amount of 
contributions (W n ) originating from all uj n 's, n = 1, . . . , oo. This is how an infinite amount of 
overtones comes into play. These overtones are multiples not of precession rate to but of the "base 
frequency" uj\ = tun / \2K{k 2 )) which is lower than cu . Hence the stresses and strains contain not 
only Fourier components oscillating at frequencies higher than the precession rate, but also components 
oscillating at frequencies lower than cu . This is a very unusual and counterintuitive phenomenon. 

The above series (called "nome expansions") typically converge very quickly, for q 1 . Note, 
however, that q — > 1 at the separatrix. Indeed, on approach to the separatrix we have: k — > 1 , 
wherefrom K{k 2 ) — > oo ; therefore q — > 1 and uj n — >• (see eqn. (|6.24|)). The period of rotation 



(see ( |6.20| )) becomes infinite. (This is the reason why near- separatrix states can mimic the principal 
one.) 

Our previous work, Efroimsky (2000), addressed relaxation in the vicinity of poles. This case 
corresponds to k 1 . For this reason we used, instead of ( |5.21| - |5.23| ), trivial approximations 



uj\ ~ uj , sn(uit, k 2 ) ~ sin(co>£) , cn(ujt, k 2 ) w cos(cjt) , dn(ujt, k 2 ) rs 1. These approximations, 
along with (|6T8l- 16.18"| ) enabled us to assume that the terms o 2 z and a 2 z in d5li| ) are associated with the 
principal frequency uj , while < o 2 x >, < o 2 y >, < a 2 zz > , < (Tra) 2 > and o xy are associated 
with the second harmonic 2 to . No harmonics higher than second appeared in that case. However, if 
we move away from the poles, parameter k will no longer be small (and will be approaching unity as 
we approach the separatrix). Hence we shall have to take into account all terms in ( |6.21| - 16.23| ) and, as 



a result, shall get an infinite amount of contributions from all oj n 's in ( |4-.6| - P^6] ). Thus we see that the 
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problem is very highly nonlinear. It is nonlinear even though the properties of the material are assumed 
linear (strains e are linear functions of stresses a ). Retrospectively, the nonlinearity originates be- 
cause the dissipation rate (and, therefore, the relaxation rate) is proportional to the averaged (over the 
cycle) elastic energy stored in the body experiencing precession-caused alternating deformations. The 
average elastic energy is proportional to < a e > , i.e., to < a 2 > . The stresses are proportional to the 
components of the acceleration, that are quadratic in the components of the angular velocity (|5.1|- |6.2j). 
All in all, the relaxation rate is a quartic form of the angular-velocity components that are expressed by 
the elliptic functions Qo\2T| - |Q3D . 

A remarkable fact about this nonlinearity is that it produces oscillations of stresses and strains not 
only at frequencies higher than the precession frequency cu but also at frequencies lower than cu. This 
is evident from formula ( |6.24| ): the closer we get to the separatrix (i.e., the closer k 2 gets to unity), the 
smaller the factor n /(2K) , and the more lower-than-co> frequencies emerge. 

A quantitative study of near-separatrix wobble will imply attributing extra factors of uj n /Q{uj n ) to 
each term of the series ( J3.6[ ) and investigating the behaviour of the resulting series (|4.1|). This study will 
become the topic of our next paper. Nevertheless, some qualitative judgement about the near-separatrix 
behaviour can be made even at this point. 



For the calculation of the dissipation rate (pTTD, the value of the average elastic energy < W > given 
by the sum (KR) is of no use (unless each of its terms is multiplied by u> n /Q(u n ) and plugged into 
Q4.1D). For this reason, the values of the terms < cr?- > entering ( |4^ ) are of no practical value either; 
only their expansions obtained by plugging ( |6.21| - 6.23[ ) into ( |5T%| - |6.18| ) do matter. Nonetheless, let 



us evaluate < W > near the separatrix. To that end, one has to calculate all < of > 's by evaluating 



( |5TT5| - I5TTED - Direct integration in ( |5TT5| - |5TT9p leads to: 



3k A 



2E 



k ) - 3 



(6.25) 
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* f 2 (* 2 - + 
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K 
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(6.26) 



iP {# (i + * 2 ) + (* 2 - 0} 



(6.27) 



(6.28) 



K and E being abbreviations for the complete elliptic integrals of the 1st and 2nd kind: 

rir/2 



K = K(k 2 ) = f (1 - k 2 sin 2 ^)" 172 di) 
Jo 

E = E(k 2 ) = f n/2 (1 - k 2 sin 2 ^) 1/2 d^j . 
Jo 



(6.29) 



> 



In the limit of k — > 1 , the expression for K will diverge and all Sj will vanish. Then all < cr 
will also become nil, and so will < W > . As all the inputs < W(to n ) > in Q4-.1D are nonnegative, each 
of them will vanish too. Hence the relaxation slows down near the separatrix. Moreover, it appears to 
completely halt on it. How trustworthy is this conclusion? On the one hand, it might have been guessed 
simply from looking at Q6.2(J| ): since for k — > 1 the period AK{k 2 ) diverges (or, stated differently, 
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since the frequencies u n in ( |6.24| ) approach zero for each fixed n, then all the averages may vanish). 



On the other hand, though, the divergence of the period undermines the entire averaging procedure: for 
r — > oo , expression Q3.3D becomes pointless. Let us have a look at the expressions for the angular- 
velocity components near the separatrix. According to (Abramovits & Stegun 1965), these expressions 
may be expanded into series over small parameter (1 — k 2 ) : 

Qi — 7 dn (ait, k 2 ^j = 7 {sech(ujt) + 
+ - (1 - k 2 ) [sinh(ut) cosh(ut) + out] sech(ut) tanh(ut)^ + O ((1 - /c 2 ) 2 ) , (6.30) 

VL2 = (3 sn (ujt, k 2 ^j = (3 {tanh (cut) + 
+ 1 (1 - k 2 ) \sinh(uot) cosh(uot) - ut) sech 2 (uot)^ + O ((1 - A; 2 ) 2 ) , (6.31) 

Q3 = a cn (cut, k 2 ^j = a {sech{ujt) — (6.32) 

- - (1 - k 2 ) [smh(uot) cosh(ujt) - cut] sech{uot) tanh{ut)\ + O ((1 - A; 2 ) 2 ) . (6.33) 

These expansions will remain valid for small k 2 up to the point k 2 = 1 , inclusively. It doesn't mean, 
however, that in these expansions we may take the limit of t — > 00 . (This difficulty arises because 
this limit is not necessarily interchangeable with the infinite sum of terms in the above expansions.) 
Fortunately, though, for k 2 = 1 , the limit expressions 

Q 1 = 7 dn(ujt, 1) = 7 sech{ujt) , (6.34) 
VL 2 = 13 m(ut, 1) = (3 tanh(ujt) , (6.35) 
Q 3 = a cn{ojt, 1) = a sechiut) (6.36) 



make an exact solution to fl2.ip. Thence we can see what happens to vector fl when its tip is right on 
the separatrix. If there were no inelastic dissipation, the tip of vector Ct would be slowing down while 
moving along the separatrix, and will come to halt at one of the middle-inertia homoclinic unstable poles 
(though it would formally take Ct an infinite time to get there, because fix and f2 3 will be approaching 
zero as ~ exp(— cut) ). When ft gets sufficiently close to the homoclinic point, the precession will 
slow down so that an observer would get an impression that the body is in a simple-rotation state. In 
reality, some tiny dissipation will still be present even for very slowly evolving Ct . It will be present 
because this slow evolution will cause slow changes in the stresses and strains. The dissipation will 
result in a further decrease of the kinetic energy, that will lead to a change in the value of k 2 (which 
is a function of energy; see ( |6.3| ) and (|6.4[)). A deviation of k 2 away from unity will imply a shift 
of Cl away from the separatrix towards pole C. So, the separatrix eventually will be crossed, and the 
near- separatrix slowing-down does NOT mean a complete halt. 

This phenomenon of near- separatrix slowing-down (that we shall call lingering effect) is not new. 
In a slightly different context, it was mentioned by Chernous'ko (1968) who investigated free precession 
of a tank filled with viscous liquid and proved that, despite the apparent trap, the separatrix is crossed 
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within a finite time. Recently, the capability of near-intermediate-axis rotational states to mimic simple 
rotation was pointed out by Samarasinha, Mueller & Belton (1999) with regard to comet Hale-Bopp. 

We, thus, see that the near-separatrix dissipational dynamics is very subtle, from the mathematical 
viewpoint. On the one hand, more of the higher overtones of the base frequency will become relevant 
(though the base frequency itself will become lower, approaching zero as the angular-velocity vector 
approaches the separatrix). On the other hand, the separartrix will act as a (temporary) trap, and the 
duration of this lingering is yet to be estimated. 

One should, though, always keep in mind that a relatively weak push can help the spinning body to 
cross the separatrix trap. So, for many rotators (at least, for the smallest ones, like cosmic-dust grains) 
the observational reality near separatrix will be defined not so much by the mathematical sophistries 
but rather by high-order physical effects: the solar wind, magnetic field effects, etc... In the case of a 
macroscopic rotator, a faint tidal interaction or a collision with a smaller body may help to cross the 
separatrix. 



VII. APPLICATION TO ASTEROIDS AND COMETS 

Let us begin with 4179 Toutatis. This is an S-type asteroid analogous to stony irons or ordinary 
chondrites, so the solid-rock value of p Q suggested in Efroimsky & Lazarian (2000) may be applicable 
to it: pQ ps 1.5 x 10 13 dyne/cm 2 = 1.5 x 10 12 Pa . Its density may be roughly estimated as 
p = 2.5 x 10 3 kg/m 3 (Scheeres et al. 1998). Just as in the case of ( |5.32| ), let us measure the time At 
in years, the revolution period T in hours (T( feonr . s ) = 175), the maximal half -radius a in kilometers 
( a {km) — 2.2 ), and 6 in angular degrees (| A8°\ = 0.01). Then (j6^) will yield: 



T 3 

At {years) ps 5.1 x iQ-2 {h ° urs) = 5.6 x 10 4 years (7.1) 

a {km) 

Presently, the angular- velocity vector Ct of Toutatis is at the stage of precession about A (see Fig. 2). 
However its motion does not obey the restriction (cos 2 9) < 1/7 under which (j6.7|) works well. A 
laborious calculation based on equations (2.16) and (A4) from Efroimsky (2000) and on formulae (1), 
(2) and (11) from Scheeres et al (1998) shows that in the case of Toutatis (cos 2 6) ps 2/7 . Since the 
violation is not that bad, one may still use ( |7.1[ ) as the zeroth approximation. Even if it is a two or three 
order of magnitude overestimate, we still see that the chances for experimental observation of Toutatis' 
relaxation are slim. 

This does not mean, though, that one would not be able to observe asteroid relaxation at all. The 
relaxation rate is sensitive to the parameters of the body (size and density) and to its mechanical prop- 
erties (pQ), but the precession period is certainly the decisive factor. Suppose that some asteroid is 
loosely-connected ( pQ = 5 x 10 12 dyne/cm 2 = 5 x 10 11 Pa and p = 2 x 10 3 kg/m 3 ), has 
a maximal half-size 17 km, and is precessing with a period of 30 hours, and is not too close to the 
separatrix. Then an optical resolution of | A0°| = 0.01 degrees will lead to the following time interval 
during which a 0.01° change of the precession-cone half-angle will be measurable: 

At (years) ~ 2.12 x 10" 2 {h ° urs) = 2 years (7.2) 

a {km) 

which looks most encouraging. In real life, though, it may be hard to observe precession relaxation of 
an asteroid, for one simple reason: too few of them are in the states when the relaxation rate is fast 
enough. Since the relaxation rate is much faster than believed previously, most excited rotators have 
already relaxed towards their principle states and are describing very narrow residual cones, too narrow 
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to observe. The rare exceptions are asteroids caught in the near-separatrix "trap". These are mimicing 
the principal state. 

On these grounds, it is easy to guess the rotational state of 433 Eros: since it is not in a sweep-tumble 
mode, then most probably it is not precessing at all, or keeps an extremely narrow residual cone. An 
almost circular precession with a half-angle of several degrees is very improbable because most likely 
it has already been transcended. Indeed, the observations have indicated no visible wobble (Yeomans 
et al 2000). 

What about comets? According to Peale and Lissauer (1989), for Halley's comet a « 
10 10 dyne /cm 2 = 10 9 Pa while Q < 100, like for the regular ice. We are unsure if the values of 
order 100 for Q are acceptable; we would be more comfortable with values close to those of firn 
(heavy coarse-grained snow): Q« 1, TherQ uQ m 10 10 dyne /cm 2 = 10 9 Pa. As for the density of 
the cometary material, it is probable that the average density of a comet does not deviate much from 
1.5 x 10 3 kg/m 3 . Indeed, on the one hand, the major part of the material may have density close to that 
of firn, but on the other hand a typical comet will carry a lot of crust and dust on and inside itself. Now, 
consider a comet of a maximal half-size 7.5 km (like that of Halley comet (Houpis and Gombosi 1986)) 
precessing with a period of 3.7 days ~ 89 hours (just as Halley doesQ). If we once again assume the 
angular resolution of the spacecraft-based equipment to be | A9°\ = 0.01, it will lead us to the following 
damping time: 

T 3 

At {years) « 5.65 x 10~ 5 {h ° urs) = 0.7 year . (7.3) 

a (km) 

This means that the cometary-relaxation damping may be measurable. 

It also follows from ( [73| ) that, to maintain the observed tumbling state of the Comet P/Halley, its jet 
activity should be sufficiently highf 9 ] 

VIII. APPLICATION TO ASTEROID 433 EROS IN LIGHT OF RECENT OBSERVATIONS 

As already mentioned in the above section, asteroid 433 Eros is in a spin state that is either principal 
one or very close to it. This differs from the scenario studied in (Black et al 1999). According to that 
scenario, an almost prolate body would be spending most part of its history wobbling about the minimal- 
inertia axis. Such a scenario was suggested because the gap between the separatrices embracing pole 
C on Fig. 2 is very narrow, for an almost prolate top, and therefore, a very weak tidal interaction or 
impact would push the asteroid's angular velocity vector Ct across the separatrix, away from pole C. 
This scenario becomes even more viable due to the "lingering effect" described in section V, i.e., due to 
the relative slowing down of the relaxation in the closemost vicinity of the separatrix. 



'Our estimate of Q still remains rough, because the inner layers of the comet may contain amorphous water 
frost (Prialnik 1999), material whose attenuation properties may differ from those of firn. 

18 Belton et al 1991, Samarasinha and A'Hearn 1991, Peale 1992 

19 The effect of outgassing upon the rotational state has been addressed in several articles. Wilhelm (1987) for 
the first time demonstrated numerically that spin states can undergo significant changes due to outgassing torques. 
This was followed by Julian (1990). A detailed numerical treatment covering effects of outgassing over many 
orbits is presented in Samarasinha and Belton (1995). 
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Nevertheless, this scenario has not been followed by Eros. This could have happened for one of the 
following reasons: either the dissipation rate in the asteroid is high enough to make Eros well relaxed 
after the recentmost disruption, or the asteroid simply has not experienced impacts or tidal interactions 
since times immemorial (since the early days of the Solar System, if we use the estimates by Burns & 
Safronov (1973) who argued that the characteristic times of asteroid relaxation may be of order hundred 
of millions to billion years). 

The latter option is very unlikely: currently Eros is at the stage of leaving the main belt; it comes 
inside the orbit of Mars and approaches that of the Earth. It is then probable that Eros during its recent 
history was disturbed by the tidal forces that drove it out of the principal spin state. 

Hence we have to prefer the former option, option that complies with our theory of precession 
relaxation. The fact that presently Eros is within less than 0. 1 degree from its principal spin state means 
that the precession relaxation process is a very fast process, much faster than believed previously^. 



IX. UNRESOLVED ISSUES 

Our approach to calculation of the relaxation rate is not without its disadvantages. Some of these are 
of mostly aesthetic nature, but at least one is quite alarming. 

As was emphasised in the end of Section II, our theory is adiabatic, in that it assumes the presence of 
two different time scales or, stated differently, the superposition of two motions: slow and fast. Namely, 
we assumed that the relaxation rate is much slower than the body-frame-related precession rate uj (see 



formulae (£23)) and (|Zq)). This enabled us to conveniently substitute the dissipation rate by its average 
over a precession cycle. The adiabatic assertion is not necessarily fulfilled when uj itself becomes small. 
This happens, for example, when the dynamical oblateness of an oblate ( J 3 > J 2 = I\ = /) body is 
approaching zero: 

(h - 1) - , h = h/I . (9.1) 

Since in the oblate case uj is proportional to the oblateness (see (5.4)), it too will approach zero, making 
our adiabatic calculation inapplicable. This is the reason why one cannot and shouldn't compare our 
results, in the limit of (h — 1) — > , with the results obtained by Peale (1973) for an almost-spherical 
oblate body. 

In the general, triaxial case, our result, should not be compared, in the limit of weak triaxiality, to 
those presented in Peale, Cassen & Reynolds (1979) and Yoder (1982), because those papers addressed 
not free dissipation but tidal dissipation. Our results, in the limit of weak triaxiality, should not be 
compared either to those obtained by Yoder & Ward (1979) for Venusian wobble-damping rate. The 
results of Yoder & Ward (1979) are correct in the limit they were designed for, i.e., for an almost 
spherical planet. None of the asteroids and comets are almost spherical; hence they are subject to our 
approach, not to that of Yoder & Ward. 

Another minor issue, that has a lot of mathematics in it but hardly bears any physical significance, is 
our polynomial approximation ( |5.13| - 15.16| , |6]8| - 16.14| ) to the stress tensor. As explained in Section V, 



this approximation keeps the symmetry Oij = cr, j and exactly satisfies ( |4.4| ) with ( |5.12| ) plugged in. The 
boundary conditions are fulfilled exactly for the diagonal components of the tensor and approximately 
for the off-diagonal elements. In the calculation of the relaxation rate, this approximation will result in 
some numerical factor, and it is highly improbable that this factor differs much from unity. 



20 Note that the complete (or almost complete) relaxation of Eros cannot be put down to the low values of the 
quality factor of a rubble pile, because this time we are dealing with a rigid monolith (Yeomans et al. 2000). 
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A more serious difficulty of our theory is that it cannot, without further refinement, give a reasonable 
estimate for the duration of the near-separatrix slowing-down mentioned in the end of Section VI. On 
the one hand, many (formally, infinitely many) overtones of the base frequency uj\ come into play near 
the separatrix; on the other hand, the base frequency approaches zero. Thence, it will take some extra 
work to account for the dissipation associated with the stresses oscillating at uj\ and with its lowest 
overtones. (The dissipation due to the stresses at these low frequency cannot be averaged over their 
periods.) 

There exists, however, one more, primary difficulty of our theory. Even though our calculation pre- 
dicts a much faster relaxation rate than believed previously, it still may fail to account for the observed 
relaxation which seems to be even faster than we expect. This paper was already in press when Andrew 
Cheng confirmed the preliminary conclusion of the NEAR team, that the upper limit on non-principal 
axis rotation is better than 0. 1 angular degree^. How to interpret such a tough observational limit on 
Eros' residual precession-cone width? Our theory does predict very swift relaxation, but it also shows 
that the relaxation slows down near the separatrix and, especially, in the closemost vicinity of poles A 
and C . Having arrived to the close vicinity of pole C , the angular- velocity vector Q must exponen- 
tially slow down its further approach to C (see the paragraph after equation (5.23)). For this reason, 
a body that is monolithic (so that its fiQ is not too low) and whose motion is sometimes influenced 
by tidal or other interactions, must demonstrate to us at least some narrow residual precession cone. 
As already mentioned, for the past million or several millions of years Eros has been at the stage of 
leaving the main belt. It comes inside the Mars orbit and approaches the Earth. It is possible that Eros 
experienced a tidal interaction within the said period of its history. Nevertheless it is presently in or ex- 
tremely close to its principal spin state. The abscence of a visible residual precession not only disproves 
the old theory but also indicates that our new-born theory, too, may be incomplete. In particular, our 
Q-factor-based empirical description of attenuation should become the fair target for criticisms, because 
it ignores several important physical effects. 

One such effect is material fatigue. It shows itself whenever a rigid material is subject to repetitive 
load. In the case of a wobbling asteroid or comet, the stresses are tiny, but the amount of repetitive 
cycles, accumulated over years, is huge. At each cycle, the picture of emerging stresses is virtually the 
same. Moreover, beside the periodic stresses, there exists a constant component of stress. This may lead 
to creation of "weak points" in the material, points that eventually give birth to cracks or other defects. 
This may also lead to creep, even in very rigid materials. The creep will absorb some of the excessive 
energy associated with precession and will slightly alter the shape of the body. The alteration will be 
such that the spin state becomes closer to the one of minimal energy. It will be achieved through the 
slight change in the direction of the principal axes in the body. If this shape alteration is due to the 
emergence of a considerable crack or displacement, then the subsequent damping of precession will be 
performed by a finite step, not gradually. 

Another potentially relevant phenomenon is the effect that a periodic forcing (such as the solar 
gravity gradient) would have on the evolution and relaxation of the precession dynamics. It is possible 
that this sort of forcing could influence the precessional dynamics of the body0. 



Andrew Cheng, personal communication. 

I am thankful to Daniel Scheeres who drew my attention to this effect. 
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X. RUBBLE HEAP VERSUS MONOLITH 



Above we mentioned one of the most important discoveries of the NEAR-Schoemaker mission: Eros 
is a well-connected monolith. This brings up an interesting issue that is still unresolved. 

At present, most astronomers lean toward the rubble-pile hypothesis, in regard to both asteroids and 
comets. The hypothesis originated in mid-sixties (Opik 1966) and became a dominating theory in the 
end of the past century (Burns 1975; Asphaug & Benz 1994; Harris 1996; Asphaug & Benz 1996; 
Bottke & Melosh 1996a,b; Richardson, Bottke & Love 1998; Bottke 1998, Bottke, Richardson & Love 
1998; Bottke, Richardson, Michel & Love 1999, Pravec & Harris 2000). 

Sometimes comets get rent apart by the tidal forces (Asphaug & Benz 1996, Sekanina 1982, Melosh 
& Schenk 1993). On these and other grounds many researchers conclude that all comets are weakly 
connected. A possible counter argument may be the following: since the comets, when warmed up 
by the Sun, are prone to tidal disintegration, then perhaps, the weakest comets have already perished 
and only the strongest have survived. Hopefully, our understanding of the subject will improve after 
the Deep Impact mission reaches its goal. Meanwhile, we would lean towards the moderate viewpoint 
(Efroimsky & Lazarian 2000): at least some of the comets are loosely connected conglomerates, but we 
do not know if all or even if most of them are like that 

In the case of asteroids, it may be unwise of us to completely reject the rubble-pile hypothesis. This 
hypothesis rests on several strong arguments the main of which is this: the large fast-rotating asteroids 
are near the rotational breakup limit for aggregates with no tensile strength. Still, we would object to 
two of the arguments often used in support this theory. One such dubious argument is the low density 
of asteroid 253 Mathilde. The low density of Mathilde (Veverka et al 1998, Yeomans et al 1998) may 
indeed evidence of high porosity. However, in our opinion, the word "porous" is not necessarily a 
synonim to "rubble-pile", even though in the astronomical community they are often used as synonims. 
In fact, a material may have high porosity and, at the same time, be rigid. 

Another popular argument, that we would contest, is the one about crator shapes. Many colleagues 
believe that a rigid body would be shattered into smitherines by collisions; therefrom they infer that the 
asteroids must be soft, i.e., rubble. In our opinion, though, a rigid but highly porous material may stand 
very energetic collisions without being destroyed, because its porous structure damps the impact. 

Finally, it is know from the construction engineering that some materials, initially friable, become 
relatively rigid after being heated up (like, for example, asphalt). They remain porous and may be prone 
to creep, but they are, nevertheless, sufficently rigid and well connected. 

For these three reasons, we expressed in Efroimsky & Lazarian (2000) our conservative opinion on 
the subject: at least some asteroids are well-connected solid chunks, though we are uncertain whether 
this is true for all asteroids. This opinion met a cold reaction from the community. However, it is 
supported by the recentmost findings. The monolithic nature of Eros is the most important of these. 
Other include 1998KY26 studied in 1999 by Steven Ostro and his team: from the radar and optical 
observations, the team inferred that this body, as well as several other objects, is monolithic (Ostro et al 
1999). 

Still, we have to admit that the main argument in favour of rubble-pile hypothesis (the absence of 
large fast rotators) remains valid. 

XL CONCLUSIONS 

1. In many spin states, dissipation at frequencies different from the precession frequency makes 
a major input into the inelastic -relaxation process. These frequencies are overtones of some "basic" 
frequency, that is LOWER than the precession frequency. Thereby we encounter a very unusual example 
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of nonlinearity: the principal frequency (precession rate) gives birth not only to higher frequencies but 
also to lower frequencies. 

2. Distribution of stresses and strains over the volume of a precessing body is such that a major 
share of inelastic dissipation is taking place deep inside the body, not in its shallow regions, as thought 
previously. These and other reasons make inelastic relaxation far more effective than believed hitherto. 

3. However, if the rotation states that are close to the separatrix on Fig. 2, the lingering effect takes 
place: both precession and precession-damping processes slow down. Such states (especially those 
close to the homoclinic point) may mimic the principal rotation state. 

4. A finite resolution of radar-generated images puts a limit on our ability of recognising whether 
an object is precessing or not. Relaxation-caused changes of the precession-cone half-angle may be 
observed. Our estimates show that the modern spacecraft-based instruments are well fit for observations 
of the asteroid and cometary wobble relaxation. In many rotation states, relaxation may be registered 
within relatively short periods of time (about a year). 

5. Measurements of the damping rate will provide us with valuable information on attenuation in 
small bodies, as well as on their recent histories of impacts and tidal interactions 

6. Since inelastic relaxation is far more effective than presumed earlier, the number of asteroids 
expected to wobble with a recognisable half-angle of the precession cone must be lower than expected. 
(We mean the predictions suggested in (Harris 1994).) Besides, some of the small bodies may be in the 
near-separatrix states: due to the afore mentioned lingering effect, these rotators may be "pretending" 
to be in a simple rotation state. 

7. Even though our theory predicts a much higher relaxation rate than believed previously, this 
high rate may still be not high enough to match the experimentally available data. In the closemost 
vicinity of the principal spin state the relaxation rate must decrease and the rotator must demonstrate the 
"exponentially-slow finish". Asteroid 433 Eros is a consolidated rotator whose Q-factor should not be 
too low. It is possible that this asteroid was disturbed sometimes in its recent history by the tidal forces. 
Nevertheless, it shows no visible residual precession. Hence, there may be a possibility that we shall 
have to seek even more effective mechanisms of relaxation. One such mechanism may be creep-caused 
deformation leading to a subsequent change of the position of the principal axes in the body. 

XII. WHAT IS TO BE DONE. 

Our further advance in the theoretical analysis of the phenomenon and in planning the appropriate 
missions should include several steps. 

1 . Our previous work (Efroimsky 2000) accounts for the dynamics at the stage when the angular- 
velocity vector f2 and the major-inertia axis of the body describe almost circular cones about the 
angular-momentum vector; that corresponds to O describing almost circular trajectories on Fig. 2. The 
next step would be to get an expression for the damping rate of a wobbling triaxial rotator at the other 
stages of precession. In particular, it would be important to estimate the duration of the near-separatrix 
lingering, i.e., the time during which a rotator can mimic a simple rotation state. 

2. Second, it is important to improve the precision of our calculation by taking into account the 
real shapes of precessing bodies: it would be more natural to model a body not by a prism but by an 
ellipsoid. This will demand a more refined mathematical approach to the appropriate boundary-value 
problem for stresses. 

3. Instrumentation on spacecraft have angular resolution of 0.01 degree (0.6 arcmin) or even better. 
It is a separate data-handling problem to make this resolution translate into a similar resolution in the 
precession-cone half-angle 9. (Or, in the general case of a triaxial or prolate rotator, into a similar 
resolution in the averaged-over-a-cycle sin 2 9 .) 
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4. Last, and by no means least, the mechanical and physical properties of the asteroid and cometary 
materials must be studied. Much work in this direction has already been done (Klinger et al. 1996; 
Muinonen & Lagerros 1998; Remo 1994, 1999; Prialnik 1999), but our knowledge of attenuation in 
small bodies still remains very basic, and consists more of hypotheses than of facts. So we are in a bad 
need of both experimental and theoretical results on attenuation in the materials asteroids and comets 
are made of. At the same time, it is the future measurements of the relaxation rate that will shed light 
on the material properties of the tumbling objects. 

5. The above program, if carried out, will open up realistic perspectives for measuring the wobbling 
bodies' relaxation rate. Observation of spin states is naturally a part of any rendezvous mission. In 
future, though, it would be better to perform not one but a series of such measurements, by each such 
mission. In the case of comets, it would be good to measure the spin state shortly before the perihelion 
(about 3 AU from the Sun or farther, i.e., before the outgassing of water begins). The second measure- 
ment should be performed shortly after the perihelion. Finally, at least one more observation would be 
in order several months (or, even better, years) later. Such a scheme would show the dynamics of both 
excitation and damping. Calculations show that wobble-damping measurements have a good chance 
of success. In the case of an asteroid, a success of such an experiment will crucially depend upon the 
structure of object: it may be difficult to observe damping of a solid-rock asteroid, because the dissipa- 
tion in solid rock is slow. In rubble-pile asteroids dissipation is several orders faster, and we may have 
a good chance of observing relaxation of such rotators. In the case of a comet, we have a very good 
chance to register precession relaxation within a less than year period of time, if the spin state is not on 
the separatrix. 
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